Phase Diagram and Incommensurate Phases in Undoped Manganites 
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We study the existence of incommensurate phases in the phase diagram of the two orbital double 
exchange model coupled with Jahn- Teller phonons and with superexchange interactions. In agree- 
ment with experimental results, we find that undoped manganites RMnOz (R being some rare earth 
element) show temperature induced commensurate-incommensurate phase transitions. In the in- 
commensurate phase the magnetic wave vector varies with temperature. The incommensurate phase 
arises from the competition between the short range antiferromagnetic superexchange interaction 
and the long range ferromagnetic double exchange interaction. 



PACS numbers: 75.47.Gk,75.10.-b. 75.30Kz, 75.50.Ee. 



I. INTRODUCTION 

Perovskites of manganese of formula (Ri_ a: A :c )Mn03 
where R denotes rare earth ions(R= La, Pr, Nd ...) and 
A is a divalent alkaline ion (A=Ca,Sr...) have attracted 
a great interest because they show a remarkable colossal 
magnetoresistance effect at doping x near one third|l|,|2(. 
In these oxides x corresponds to the concentration of 
holes moving in the e g orbital bands of the Mn ions that 
ideally form a cubic structure. From the basic point of 
view these materials are a challenge for both theoreti- 
cians and experimentalists as they show a very rich phase 
diagram 3] . As function of temperature and hole dop- 
ing, these systems present orbital, charge or spin order, 
and in a large portion of the pha se diagram these orders 
coexist H HB ES II ES IH III El El III Nanophase 
separation near a;=l/3[lg and commensurate incommen- 
surate transition near half doping seem to occur in colos- 
sal magnetoresistance manganites [TtI fljl Il9ll20| . 

Many properties of manganites depend on the competi- 
tion between the kinetic energy tending to delocalize the 
carriers and localization effects such as the Jahn- Teller 
(JT) coupling and the antiferromagnetic (AFM) coupling 
between the Mn core spins. Therefore the properties 
of manganites at intermediate doping can be described 
within a band structure picture, where the itinerant e g 
carriers have a strong ferromagnetic interaction with the 
core tig Mn spins, and are coupled with the Jahn- Teller 
distortions of the oxygen octahedra surrounding the Mn 
ions [3l ITU . However, the parent compounds, RM11O3, are 
always insulator and their physical properties were typi- 
cally described in the picture of stron g l y co rrelated Mott 
localized d-electrons |Hl 111 El 111 M 113 • However 
some spin and orbital ordering found experimentally in 
undoped materials 28, 29, 30] can not be easily described 
in this picture of strongly localized electrons. 

In ref.jH the authors examine the magnetic and or- 
bital order in a series of RMnC>3 as a function of the 
ionic radius (tr) of the rare earth ion R. For small ionic 
radius the manganites have a antiferromagnetic spin or- 



der of type A coexisting with a (jr, tt, 0) orbital ordering, 
whereas for larger values of tr the magnetic order is of 
type E. In the A phase, a Mn spin is ferromagnetically 
coupled with the Mn spins located in the same plane 
(x — y), and antiferromagnetically with the Mn spins be- 
longing to different planes. In the E phase the x — y lay- 
ers are antiferromagnetically coupled, but the magnetic 
structure within the planes is that of ferromagnetic zigzag 
chains coupled antiferromagnetically. The horizontal (x) 
and vertical (y) steps of the zigzag chains contain two Mn 
ions. For values of rji close to the critical value where 
the magnetic order changes from A-type to -E-type, the 
manganites develop different magnetic incommensurate 
phases when increasing temperature. 

It has been suggested recently [^H El E3- that the 
complete nesting between the two e g bands that occurs 
in the A structure produces a spin-orbital ordering and 
opens a gap in the energy spectrum of undoped RMnC>3. 
Quoting ref.[33j, we do not claim that the real RM11O3 
systems can be fully described by a weak coupling ap- 
proach as correlation effects can be important, although 
a treatment based in band structure calculation may be 
very useful to understand some properties of these mate- 
rials. In particular, in ref.|U, using a two orbital double 
exchange model, it was obtained that the experimental 
observed i?-phase exits in a wide region of parameter 
space, and it is adjacent to the A-type phase. 

One of the issues that remains to be understood is 
the microscopic origin of the commensurate phases ap- 
pearing near the A-type to .E-type magnetic transition. 
The aim of this work is to explain these phases using a 
realistic microscopic model. The Hamiltonian we study 
describes electrons moving in two e g bands, that are fer- 
romagnetically strongly coupled to the Mn core spins as 
well to the Jahn- Teller phonons. In addition, we also 
consider a direct superexchange interaction between the 
core Mn spins. Starting from this Hamiltonian we de- 
rive a functional that describes a temperature induced 
commensurate-incommensurate transition similar to that 
observed experimentally. 
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The main result of this work is that near the A to 
E phase transition, the competition between the near- 
est neighbor antiferromagnetic superexchange interaction 
and the double exchange induced long range ferromag- 
netic interaction, results in the appearance of incommen- 
surate phases. These phases consist of a periodic array 
of domain walls. 

The rest of the paper is organized as follows. In Sec. II 
we describe the microscopic model and we present the 
zero temperature phase diagram. In Sec. Ill we outline 
the method for obtaining the critical temperatures and 
we present the phase diagram composed of the different 
uniform phases. In Sec. IV we develop the functional for 
describing spatially modulated phases. Also in Sec. IV we 
study how the phase diagram of manganites at x—1 is al- 
tered when soliton incommensurate phases are taken into 
account. We finish in Section V with a brief summary. 

II. MICROSCOPIC HAMILTONIAN AND ZERO 
TEMPERATURE PHASE DIAGRAM. 
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FIG. 1: Zero temperature phase diagram for x=l for the 
two-dimensional DE two orbital model with cooperative Jahn 
Teller phonons. The symbols OO and OD stand for orbital 
ordered and disordered respectively. A, E and G name the 
different magnetic order defined in the text. 



We are interested in the transition between the A and 
E phases. In these phases the x — y planes are cou- 
pled antiferromagnetically and therefore we can analyze 
the properties of these states and the transition between 
them by studying a Hamiltonian which describes elec- 
trons moving in the x — y plane. The Coulomb inter- 
action between electrons prevents double occupancy and 
aligns the spins of the d orbitals. The crystal field splits 
the Mn d levels into an occupied t 2g triplet and a doublet 
of e g symmetry where 1 — x electrons per Mn have to ac- 
commodate. The Hund's coupling between the spins of 
the e g electrons and each core spin is much larger than 
any other energy in the system, and each electron spin is 
forced to align locally with the core spin texture. Then 
the e g electrons can be treated as spinless particles and 
the hopping amplitude between two Mn ions is modu- 
lated by the spin reduction factor, 

fi2 = cos — cos — + e 1 ^ 1 Wsin-sm- (1) 

where {di : 4>i] are the Euler angles of the, assumed clas- 
sical, Mn core spins IS;} , This is the so called double 
exchange (DE) model|3ll35ll3r|. 

We study a double exchange model coupled to Jahn- 
Teller (JT) phonons. We also include the antiferromag- 
netic coupling between the Mn core spins J af- 

i,j,a,a f 

+ j af SiSj + lT,(PQii + QM + Ql) 

<i,j> i 

iPi + Q2iT xi + Q 3 iT zl ) , (2) 

i 

here Cf a creates an electron in the Mn ion located at 
site i, in the e g orbital a (o=l,2 with l=|x 2 — y 2 > and 



2=\3z 2 — r 2 >). The hopping amplitude t aa , is finite for 
next neighbors Mn and depends both on the type of or- 
bital involved and on the direction u between sites i and j 
{iff = ±s/Zf x f = ±VZt x 2 f = it x 2 f = t)H. t is taken 
as the energy unit. The forth term couples the e g elec- 
trons with the three active MnC>6 octahedra distortions: 
the breathing mode Qu, and the JT modes Q21 and 
Qzi that have symmetry x 2 -y 2 and 3z 2 -r 2 respectively. 
Qu couples with the charge at site i, pi — ^2 a C^ a Ci^ a 
whereas Qu and Q^i couple with the x and z orbital pseu- 
dospin, T xi = Cf l C i2 + Cf 2 Cn and t zi = C^Cn - Cf 2 C i2 , 
respectively. The third term is the elastic energy of the 
octahedra distortions, being (3 >2 the spring constant 
ratio for breathing and JT-modes|37(. In the perovskite 
structures the oxygens are shared by neighboring MnC>6 
octahedra and the Q's distortions are not independent, 
cooperative effects being very important j^- In order to 
consider these collective effects, we consider the position 
of the oxygen atoms as the independent variables of the 
JT distortions. 

For a given value of the parameters A and Jaf, and 
a texture of core spins {S^}, we solve self-consistently 
the mean field version of Hamiltonian (01 and obtain the 
energy, the local charges {pi}, the orbital pseudospin or- 
der {r X i,T Z i} and the oxygen octahedra distortions Q a ^. 
These quantities are better described by their Fourier 
transforms, that are represented by the same symbol with 
a hat: p(G), Qi(G), ... 

In Fig.l we present the phase diagram obtained by 
solving self-consistently Eq|2| for the parent compound 
RM11O3. For the range of parameters studied, we do 
not find any solution showing charge modulation. In all 
the phases there is an electron located on each Mn ion, 
therefore in our model any gap in the energy spectrum 
is due to the spatial modulation of any other physical 
quantity. 
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For small values of Jaf the ground state is ferromag- 
netic, ^4-order. In absence of Jahn- Teller coupling this 
phase is metallic, however, for A ^0, and due to the 
perfect nesting between the e g bands, the A phase de- 
velops a gap at the Fermi energy. The Jahn- Teller cou- 
pling produces and orbital order characterized by a finite 
Fourier component of the cc-component of the pseudospin 
t x (it, tt) = t x (—tt, — 71") 7^ 0, see Fig. 2a. The orbital order 
(00) is produced by an ordered distribution of the oxy- 
gen octahedra distortions Q,2(tt,tt) = Q2(— tt, — 7r) =/= 0, 
that depends on the value of A. The amplitude of the 
distortions are modulated in order to minimize the elas- 
tic energy of the cooperative Jahn- Teller distortions, and 
the signs arise from cooperative effects. In this phase the 
(tt, tt) orbital modulation opens a gap at the Fermi en- 
ergy and the x=l manganite is an insulator, being the 
energy gap proportional to the value of the Jahn- Teller 
coupling. 

For large value of Jaf and A the system presents a G- 
type antiferromagnetic ground state and an orbital order 
characterized by a Fourier component of the pseudospin 
t x (tt, tt) = f x {— tt, —tt) 7^ 0. Each Mn ion is coupled anti- 
ferromagnetically with its next neighbors and the double 
exchange mechanism precludes the motion of the carri- 
ers, being this phase an insulator. The minimal value of 
Jaf for the occurrence of this phase depends on A, but 
in general is very large, so that this phase is rather unlike 
to occur in manganites. 

For intermediates values of Jaf the system develops 
a magnetic order of .E-type; the E phase consists of fer- 
romagnetic zigzag chains coupled antiferromagnetically. 
The horizontal and vertical steps of the chain contain two 
Mn ions. For large enough values of the Jahn- Teller cou- 
pling the E magnetic order coexists with an orbital order 
similar to the one occurring in the FM and G phases; 
this order opens a gap at the Fermi energy. In the E- 
OO phase the magnetic order is characterized by a pe- 
riodicity (2a, 2a), being a the lattice parameter of the 
square lattice. For small values of A the E phase does 
not present orbital order, although it has a gap at the 
Fermi energy. In the orbital disorder (OD) iJ-phase, the 
dispersion energy for the e g electrons along the FM zigzag 
chain is given by[U, St = (2/3)(±cos fc± Vcos 2 /c + 3), 
indicating the existence of a large band gap at occupan- 
cies corresponding to x—1. The physical origin of this 
gap is the dependence of the tunnelling probability on 
the spatial direction, t?L v = — „ for /i 7^ v. It pro- 
duces a periodicity in the hopping amplitude along the 
zigzag chain, leading to a periodic potential for the e g 
electrons. It is important to note that, contrary to the 
af=l/2 case |3 13 1 this modulation in the hopping am- 
plitude does not produce an orbital order. The E-OD 
phase is stable due to the spatial modulation of the co- 
herence between next neighbors Mn ions along the zigzag 
chain, < Cf^C i+ i,u >= - < Cf +l41 C i+2 ^ > for fi ^ v. 
The phases E-OO and E-OD have different symmetry 
and therefore the transition between this phases that oc- 
curs at finite values of A is a discontinuous transition. 
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FIG. 2: (Color online) Orbital and spin order of x=l man- 
ganites in the x-y plane. Elongated orbitals along the x (y) 
directions represent d 3x 2_ r 2 (d 3y 2_ r 2) orbitals. Circles rep- 
resent the Mn ions in a orbital disordered phase, a) Orbital 
order presents in the ferromagmetic A-OO phase, b) Same 
than a but for the E-OO phase, c) Spin order in the E-OD 
phase. The solid and the dashed lines joining the Mn ions 
indicate the modulation of the electronic coherence along the 
zigzag chains, d) Same than a but for the G-OO phase. In all 
the phases there is not modulation of the electric charge and 
there is an electron located at each Mn ion. The vectors in 
the different schemes, represent the spatial periodicity in the 
different phases, (a, a) in the A-OO and G-OO phases and 
(2a, 2a) in the E-OO and E-OD phases. In all the figures 
open and close symbols represent up and down spins. 



III. FINITE TEMPERATURE MAGNETIC 
PHASE DIAGRAM. 

The simplest way of obtaining information on the 
phase diagram corresponding to a given microscopic 
Hamiltonian is by means of the mean field approxima- 
tion. This approximation is insufficient for describing 
second order transitions, but it is successful in describing 
the phases away from the transition and in predicting 
the topology of the phase diagram. In this approach we 
compute the magnetic critical temperature of the differ- 
ent phases. 



A. Ferromagnetic (A) phase 

In this phase all the Mn spins point, on average, in a 
particular direction, and there is a finite relative magne- 
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tization <m>. Using a virtual crystal approximation, we 
consider a unique value for the spin reduction factor foj 
that corresponds to its expectation value jTRIs^. 



J -< 



1 + cos Oi 



1+ < cos On > 



l+<m> 2 



(3) 

A reduction of < m > produces a decrease in fij and 
therefore in the kinetic energy. In this way the impor- 
tance of the Jahn- Teller coupling increases as the tem- 
perature decreases. The internal energy per Mn ion of 
this phase can be written as 



E A = e A (\, <m>) + 2J AF <m> 2 



(4) 



where the electronic energy per Mn ion, sa(\ <m>), 
depends in a complicated way on A and <m> and has to 
be obtained numerically by solving Eq.J2J). 

In order to describe thermal effects it is necessary to 
compute the free energy. As the entropy of the carriers 
is very small [Tij we only include the entropy of the clas- 
sical Mn spins. We use a mean field approximation that 
neglects spatial correlations and assume for each indi- 
vidual spin a statistical distribution corresponding to an 
effective magnetic field llj,|36j. In this molecular field ap- 
proximation the entropy of the Mn spins takes the form, 



S{<m>) 



log 2 



<m>* 



20 



<m>^ 



(5) 



Using this expression for the Mn spins entropy the total 
free energy of the system for small values of <m> takes 
the form 

F{<m>) = E A -TS(<m>) 
And the Curie temperature of the A phase is, 



T c = -- 



2 de A (\,<m>) 



d <m> 2 



<frj>=0 



4 

3 JAF 



(6) 



For finite A the derivative has to be calculated numeri- 
cally. From higher derivatives of the internal energy with 
respect to the magnetization we obtain that the transi- 
tion is second order. In Eq|S|we notice that for a given 
value of the Jahn- Teller coupling the Curie temperature 
decreases linearly with the superexchange antiferromag- 
netic coupling Jaf- 



B. Antiferromagnetic E phase. 

The magnetization of the E phase is described by the 
relative amount of saturation in each zigzag chain <m s >. 
In the virtual crystal approximation fluctuations are ne- 
glected and the hopping is modulated by the spin reduc- 
tion factor that is different along the zigzag FM chain, 



f FM than between the AFM coupled chains, f AF [1113 



f FM (<m s >) = 



f AF (<m s >) = 



1+ <m s > 2 



I- <m s > z 



(7) 



The internal energy of this phase depends on A, and < 
ms>, and can be written as, 



E E = s E {\<m s >) 



(8) 



As each Mn spin core is surrounded by two Mn spins 
coupled FM and other two coupled AFM, the superex- 
change energy is zero. 

In order to compute the Neel temperature of the E 
phases, we introduce an effective field for each spin sub- 
lattice. Taking into account that both, the magnetiza- 
tion and the effective magnetic field, have different sign 
in each sublattice, we end up with the same expression 
for the entropy than in the FM phase, but just changing 
<m> by <ms> 15]. With this the free energy takes the 
form, 

F(<m s >,A) = E E -TS{<m s >) 



F(0,A)+ <m s > z U T + a 



with 



E l 



cte + a <ms> 



-b <m s > 4 +c <m s > 6 



(9) 



(10) 



Due to the symmetry of the E phase the coefficient a 
is zero, and the Neel temperature depends on the coeffi- 
cients b and c. Numerically, the coefficient b is negative 
and the transition from the E to the paramagnetic (PM) 
phase is a first order phase transition. 

It is interesting to analyze the origin of the nega- 
tive sign of the quartic term. In the OD-case the 
Jahn- Teller coupling is not large enough to produce 
orbital order. In this situation the electronic energy 
is just kinetic energy. Therefore near < ms >—0 we 
would expect that the electronic energy could be ob- 
tained perturbatively from the paramagnetic energy as, 
« ^(f FM (<m s >) + f AF (<ms>))e° E . Here 
£b(0,0) is the paramagnetic energy per Mn ion. 



Expanding the spin reduction factors near <mj> we 
find £7f =0 w (l - | <m s > 4 <m s > 8 -...) e%. As 

the electronic energy of the paramagnetic phase is nega- 
tive, the last expression suggests that the Neel temper- 
ature should be zero. However, numerically, we find a 
finite Neel temperature even for A=0. This discrepancy 
occurs because, as commented above, in the E-OD phase 
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the minimization of the kinetic energy produces a mod- 
ulation of the electron coherence along the zigzag chain, 
< Cf^C i+ljV >=- < Cf +1 ^C i+2 ,v > for \i ^ v. We de- 
scribe this modulation by a order parameter £ that rep- 
resents the (§, §) Fourier component of the electron co- 
herence. This order parameter is coupled with the stag- 
gered magnetization and the functional describing the 
electronic energy has the general form, 



E 



A=0 



1 

71 



(f FM (<™s>) + f AF (<ms>)) 



a£ 2 + ft <m s > 2 



(11) 



where we have included the elastic energy associated with 
the electron coherence and the minimal coupling between 
the staggered magnetization and the electron coherence. 
Minimizing this energy with respect the coherence pa- 
rameter £, we find £ = — < nig > 2 . Introducing this 
value in the expression of the electronic energy, Eq.JTTJ, 
we obtain, 



, 1 

1 — - <ms > 



<m s > +■ 



(12) 



and for strong enough coupling between <mj> and the 
orbital coherence, the quartic term is negative and a fi- 
nite Neel temperature is expected. It is therefore the 
coupling between the electron coherence and the stag- 
gered magnetization the responsible for the occurrence 
of a finite Neel temperature. 

In the OO-case there exits a finite orbital order param- 
eter Ta;(7r, tt), that is coupled with the staggered magne- 
tization and it is the responsible of the existence of finite 
Neel temperature. 
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FIG. 3: Phase diagrams T-Jaf for the two dimensional 
DE two orbital model with cooperative Jahn- Teller phonons 
and 2=1. In a) we plot the A = case and in b) the 
A = 1.2t case. Continuous lines represent first order transi- 
tions whereas dashed lines indicate second order transitions. 
The abbreviations naming the different phases are explained 
in the text. 



C. Temperature- Jaf magnetic phase diagram 

In Fig. 3, we plot the T- Jaf magnetic phase diagrams 
for A=0 (a) and \=\.2t (b). These phase diagrams have 
been obtained by minimizing and comparing the free en- 
ergy of the A, E and paramagnetic phases. For A=0, 
all the phases are disordered in the orbital sector, how- 
ever, for large enough values of A, the E and A phases 
present orbital order. In the later case, A=1.2t, we find 
that the critical temperature associated with the orbital 
order is much larger than the magnetic critical temper- 
atures and, therefore, in Fig. 3b the paramagnetic phase 
presents orbital order. In any case, it is important to 
note that, from the magnetic point of view, both phase 
diagrams are topologically equivalent. At large temper- 
atures the systems are always paramagnetic, for small 
Jaf and small temperature the systems present ferro- 
magnetic order, whereas for small temperature and mod- 
erates values of Jaf an antiferromagnetic order of type 
E appears. For very large values of the AFM coupling, 
not shown in Fig. 3, antiferromagnetic order of type G 
would appear. The Curie temperature corresponding to 
the paramagnetic- A phase transition decreases linearly 



with Jaf, Eq.6, until it reaches the, Jaf independent, 
Neel temperature corresponding to the paramagnetic-i? 
transition. As discussed in the previous subsection the 
A-paramagnetic transition is second order while, because 
of the coupling between different order parameters, the 
.©-paramagnetic transition is first order. 

The phase diagrams present a Lifshitz point where the 
uniform ferromagnetic A phase, the modulated ordered 
E phase and the paramagnetic disordered phase meet. 
Near the Lifshitz point there is a range of values of Jaf 
where, by increasing the temperature, the system under- 
goes an E-A transition followed by an ^4-PM transition. 
The topology of this phase diagram is similar to that 
of a Ising model with competing interactions. In that 
model, near the Lifshitz point, solitons, spatially mod- 
ulated phases and commensurate incommensurate tran- 
sitions appear when the temperatures varies [io). In the 
next section we explore the possible existence of solitons 
and incommensurate phases in the model described by 
the Hamiltonian Eq[21at x=l and near the Lifshitz point 
that appears in the T-Jaf phase diagram, Fig|3| 
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FIG. 4: (Color online) Spin order of x=l manganites in the 
x-y plane, (a) Corresponds to the E phase whereas (c) repre- 
sents the A phase. The direction perpendicular to the zigzag 
chains is shown in (a). The averaged magnetization along the 
z direction for the E and A phases are plotted in (c) and (d) 
respectively. 



IV. SOLITON THEORY AND SPATIALLY 
MODULATED PHASES. 

A. Landau functional 

Magnetically, the phases described in the previous sec- 
tions only vary along the direction perpendicular to the 
chains, see Fig0 In the mean field approximation the 
E phase is described by a spin density wave of the form 
<S>— \[2mo cos((7o-z + tt/4), with qo = ir/2 and here z 
is the position of the atoms along the direction perpen- 
dicular to the chains. We are taking the distance between 
first neighbors diagonal lines of atoms, y/2a/2, as the unit 
of length. In general the expression 



< S >— V~2 mo cos [q a z + — + 9(z] 



(13) 



describes different spatially modulated magnetic phases. 
With 9(z)—0 it describes the E phase, whereas with 
6{z)=— qqz, it represents the average magnetization in 
the position independent ferromagnetic A phase. The 
case mo=0 corresponds to the paramagnetic case. In 
general 9(z) describes phases where the average magne- 
tization changes along the direction perpendicular to the 
chains. 

Both the E and A phases are commensurate with the 
underline lattice, here we are going to study the existence 
of solitons and incommensurate phases in the system. 
The solitons are static domain walls between commen- 
surate domains. In this section we present a formalism 
which makes contact with phenomenological theories and 



provides the basis for the calculation of the nature of the 
phase diagram at finite temperature. 

We want to build a Landau theory functional where 
the order parameter is the modulation of the average spin 
along the diagonal direction. Following reference |4lll42|. 
in Ea ll3l we consider the amplitude of the magnetization, 
mo, as constant in the space. 

In order to set up a Landau theory we need to calcu- 
late the different contributions to the free energy. For 
a magnetization given by Eci ll3l and near the order dis- 
order magnetic transitions, the entropic contribution to 
the free energy can be estimated as described in Sec. Ill, 



TS~k B T 



log2 



27 
40 ? 



99 
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9 297 



dz (14) 



The superexchange antiferromagnetic interaction takes 
the form 



E A f — %Jaf ml J sin (V0(z)) dz . 



(15) 



being V9(z) the derivate of 9(z) with respect z. In the 
previous expression we have treated the position z as a 
continuous variable and we have discarded second and 
higher derivates of 9 with respect the position. 

Concerning the electronic contribution to the internal 
energy, E e , we assume that it is local an can be written 
as E e = J dz £(z), being £(z) the electronic energy density. 
We expect that it can be expanded in powers of the order 
parameter < S > and its derivates, 

£(z)=£° + ai(z) <S(z)> 1 +bi{z)(V<S(z)>) 1 

+ Q, m (z) <S(z)> 1 (V<S(z)>) m -(16) 

Here the sum over repeated indices is assumed, and be- 
cause the symmetry of the system only even powers of 
<S> and W<S(z)> contribute. The coefficients 5/, 6/ 
and ci tm are periodic in z with the periodicity of the crys- 
tal lattice and for a magnetization of the form Eq^| the 
density of electronic energy can be written as, 



£{z) 



£ +a 2 m - 
+ (64TO0 
+ (c 2 ml 



O4TO 



b 6 ml) cos 49(z) 



C4TOQ 



+ c 6 m 6 )(V9(z) + q ) 2 . (17) 



Here we have neglected higher terms in the derivates of 
the phase 9 and, as there are first order phase transitions 
in some part of the phase diagram, we keep terms up 
to the sixth power in mo. The second term in Ea ll7l is 
the umklapp term that favors the modulated commensu- 
rate solutions, 9(z)=0, -|, tt, 3-|, corresponding to the 
E phase. The last term is an elastic energy which favors 
the occurrence of the ferromagnetic A-phase, 9(z)=— qoz. 
The competition between the elastic and the umklapp 
will produce the existence of solitons and incommensu- 
rate phases. 
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^From the expression of the electronic energy of the 
i?-phase as function of the order parameter to , we no- 
tice 02=— C2<Zq- Analyzing the dependence of the elec- 
tronic energy on a constant phase (9(r) = 9q), we find 
&4=a4 + C4 qg and be=ae + cq q^. In this way only two 
subsets of parameters (for instance 6's and c's) are inde- 
pendent. Finally, for each value of A, we perform micro- 
scopic calculations of the electronic energy of the A and 
E phases and obtain the numerical values of the coeffi- 
cients 6's and c's respectively. 

Adding the entropy, Ea ll4l the antiferromagnetic en- 
ergy, Eg 1151 and the electronic internal energy, Ea ll7l we 
obtain the following expression for the free energy of the 
system, 
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F = F (T,m ) 
+ cj \ (V9{z) 



9o j 



w(l + cos46»(z)) 



dz (18) 



with, 
F (T,m ) 



, „ 3 2 9 4 198 , , , . 
\ og 2 + -ml + — ml + — ml)T 



2J AF ml +e° e - c 2 mlql 



C 



2c 2 



,Jaf\ 2 
^—JT m o 
% J 



2ce,m 







and 



w 



(64 Z M<± & - WqT) 

c 



(19) 



(20) 



(21) 



In the limit w — ► 0, the elastic contribution is the 
more important term and the phase 9{z) tends to be 
6{z) = — <7o z. On the contrary, for large values of w, the 
umklapp term is dominant and 9 wants to get a constant 
value, 9—0, 7r, 3^. A transition between the commen- 
surate phase, 9 = and the uniform ferromagnetic phase 
takes place because of the competition between these two 
terms; by tuning the values of Jaf and T we are going 
to see that a soliton incommensurate phase appears be- 
tween these two limits. 

For a given temperature and a particular value of Jaf, 
the constant amplitude mo and the phase function 9(z) 
that characterize the solution are obtained by minimizing 
the functional Eq^] For each too , the phase 9(z) should 
satisfy the sine-Gordon equation, 



FIG. 5: Phase diagram T — Jaf as obtained by minimizing 
the free energy Eq |18l The parameters of entering in the free 
energy are obtained by minimizing the microscopic Hamilto- 
nian Eq|5|for A=0 and x=l. Continuous lines represent first 
order transitions whereas dashed lines indicate second order 
transitions. The shadow region indicates the region where the 
incommensurate phase exits. 



In general the solutions of Eq^| is a soliton lattice 
formed by a regular array of domain walls, L. At each 
soliton the phase 9 tumbles \. The deviation of the av- 
erage wavevector q from go is inversely proportional to 
the distance between the domain walls, 
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7T 

2L 



(24) 



The value of q is proportional to the soliton density and 
is obtained by minimizing the free energy following the 
procedure outlined in references |40l l4ll l42l |43 | 

In the soliton lattice phase the magnetic periodicity 
along the z-direction is characterized by the wavevector 



q 2 2L 



(25) 



In the i5-phase there are not solitons in the system L=oo 
and the wavevector of the magnetic modulation is 
In the continuous approximation the ferromagnetic A- 
phase corresponds to a extremely dense lattice soliton, 
L=l. In this limit 9 is too quickly varing, the continuous 
approximation is not valid and we take the criterium that 
for L < 1.1, the soliton lattice is the ferromagnetic A- 
phase. 



ld 2 9 



4w sin A9 = 



2dz 2 

which has soliton-like solutions of the form, 



9(z) = tan exp (A^/wz) 



(22) 



(23) 



This solution is a domain wall which separates two almost 
commensurate z regions. 



B. Results 

In order to find inhomogeneous phases in manganites 
at x=l, we have minimized the free energy Ea ll8l with the 
coefficients a's, &'s and c's obtained from the microscopic 
model described in Section II. The solutions are char- 
acterized by the value of the magnetization too and the 
density of solitons q. We present results for the case A=0, 
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FIG. 6: (Color online). Temperature dependence of the mag- 
netic wave vector, q, and of the order parameter amplitude, 
mo (inset), for two values of Jaf- The wavevector q=\ cor- 
responds to the _E-phase and q=0 to the A-phase. The results 
are obtained by minimizing the free energy (Eq. I18H . 



but similar results are obtained for finite Jahn- Teller cou- 
pling. 

If we consider only the uniform solutions, E and A 
phases, the minimization of the free energy results in 
the phase diagrams already presented in Sec. Ill, FigO 
When inhomogeneous solutions are considered, we ob- 
tain the phase diagram shown in Fig0 Several com- 
ments on this phase diagram are in order, (i) There is 
not solitons for values of Jaf larger than the antiferro- 
magnetic coupling corresponding to the Lifshitz point. 
The paramagnetic- E phase transition is first order, with 
a large jump in the value of m<?, and therefore in the 
value of w. In this situation, large values of w, the umk- 
lapp term is much stronger than the elastic term and 
the system prefers to be commensurate with the lattice. 
Note than in the paradigmatic Ising model with compet- 
ing interactions 40] all the transitions are second order 
and incommensurate phases appear at both sides of the 
Lifshitz point, (ii) For small values of Jaf the elastic 
term is very strong and the solution corresponds to a 
dense soliton phase. For small values of Jaf the dis- 
tance between solitons is smaller than the cutoff and we 
consider that this commensurate phase is actually, in the 
discrete real crystal, the ferromagnetic A-phasc. (iii) For 
intermediate values of Jaf, the competition between the 
elastic and the umklapp term results in the appearance 
of incommensurate solitonic phases. 

In Fig|Sl the shadow region indicates the incommen- 
surate phase. The frontier of this phase with the ferro- 
magnetic A-phase is diffuse because, as we have already 
discussed, the difference between a dense soliton phase 
and the ferromagnetic A-phase is based in a criterium on 
the distance between solitons. Typical temperature de- 



pendence of the magnetic wavevector and magnetization 
amplitude near the incommensurate phase is illustrated 
in FigEl For low temperature the system is in the com- 
mensurate -E-phase, corresponding to a wave vector q=%. 
At low temperatures the spins are highly polarized and 
that makes the magnetic modulation too rigid to allow 
solitons. As temperature increases, the amplitude of the 
spin-wave, mo decreases and, at a temperature that de- 
pends on Jaf, a jump to the solitonic phase takes place. 
For values of Jaf closer to the Lifshitz point, the incom- 
mensurate phase appears at temperatures close to the 
A-paramagnetic critical temperature. In that case the 
amplitude of the magnetization mo is small and hence 
the magnetic wavevector of the incommensurate phase 
is also very small. Therefore in this part of the phase 
diagram the incommensurate phase is similar to the A- 
phase. 

For values of Jaf near the zero temperature A to E 
phase transition, there is a small portion of the phase 
diagram where by decreasing the temperature the sys- 
tem evolves first from a paramagnetic, q — 0, phase to 
a incommensurate phase characterized by a finite q and 
then to a ferromagnetic A phase without magnetic mod- 
ulation, q = 0. 

The magnetic phase diagram shown in Fig'0 contains 
the essence of the magnetic properties experimentally ob- 
served in undoped manganites|2^|. Systems with large 
hopping amplitude, {Jaf ft small) as LaMnC>3 have a 
ground state with a magnetic order of type A and a 
relatively large Neel temperature. The relative value of 
the AFM coupling increases when the ionic radius of the 
rare earth in the manganite increases. Therefore we un- 
derstand the experimental decrease of the Neel tempera- 
ture in the series of RMn03 (R=La, Pr, Nd, Sm) as the 
diminution of the ^-paramagnetic critical temperature 
when Jaf increases, see FigO and FigO Experimen- 
tally it is observed that for large enough ionic radius, 
HoMnC>3, the ground state of the undoped manganite 
has a magnetic order of type E, and present incommen- 
surate phases when temperature increases. We claim that 
this situation corresponds in the phase diagram FigEJto 
values JafI^ in the range 0.18-0.20. Experimentally is 
also observed that in some compounds as TbMnC>3 and 
GdMnOs, when the temperature increases the system 
undergoes two phase transitions, first a a ferromagnetic- 
incommensurate transition and at higher temperatures a 
incommensurate -paramagnetic transition. In the phase 
diagram presented in FigEJthis behavior occurs for values 
of Jaf near 0.18i. 



V. SUMMARY 

By starting from a microscopic Hamiltonian we have 
derived an expression for the free energy of undoped man- 
ganitcs. Using a realistic model we have quantified the 
competition between the short range superexchage an- 
tiferromagnetic interaction, and the long range double 
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exchange ferromagnetic interaction. The competition be- 
tween these interactions results in the existence of mag- 
netic incommensurate phases as the recently experimen- 
tally observed in undoped manganites. The incommen- 
surate phases can be described as arrays of domain walls 
separating commensurate phases by a distance that de- 
pends on temperature. The results presented in Fig[5] 
explain qualitatively the experimental results published 



in ref. 29] 
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